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Abstract. The Random Field Ising Model (RFIM) is the simplest physical model 
reflecting effect of quenched disorder on the different types of phase transitions in 
solids. The presence of multiple energy minima in the RFIM is an important feature 
determining main physical properties of the systems with a field type of quenched 
disorder. In particular, according to computer simulations of Sethna at al. performed 
for a zero temperature RFIM, the irreversible loop of magnetization appears as direct 
result of their presence. In the present paper the analytic method based on the total 
number of energy minima estimation is developed to study the hysteretic behavior of a 
zero temperature RFIM. As a result, the magnetization hysteresis loop is obtained 
as a marginal curve separating metastable and completely unstable states in the 
magnetization-external field plane. 
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1. Introduction 



The Random Field Ising Model (RFIM) is the simplest physical model reflecting the 
influence of different type of defects or quenched disorder on the different types of 
phase transitions in solids. Starting from the pioneering work of Imry and Ma [Q a 
great number of investigations have been made to understand the nature of different 
phase transitions in presence of random bond or field type of quenched disorder 
0, 1^, |, ^ H, 1^, 1^ . The literature contains several reviews on the theory and experiments 



IP, |T0|, |Tl|, [T^, |T3|, |T^, |T^. The investigations developed in time from the equilibrium 
properties and lower critical dimensionality to dynamics of domain walls, metastability 
and hysteretic properties. The RFIM is now also used to explain a role of impurities 
in displacive transitions and random fields effect in charge-density waves. The multiple 
energy minima in the RFIM is an important feature determining main equilibrium 
and non-equilibrium physical properties of these systems. The most important effect 
resulting from presence of the infinite number of local energy minima in RFIM can be 
expected out of equlibrium (for example, in zero temperature case) when the thermal 
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fluctuations become negligibly small. In particular, the irreversible magnetization loop 
appears as direct result of multiple energy minima in RFIM. This is confirmed by the 
direct computer simulations of magnetization processes reported in several publications 
(see |T^, jl^)- Due to great number of possible stationary microscopic configuration the 
magnetization process can be considered as an infinite bifurcation sequence of transitions 
between the multiple stationary (metastable) states when the external field is changed. 
Such a transitions occur each time as the external field approaches the corresponding 
critical value where the current spin configuration becomes unstable. Then, this unstable 
configuration is changed at a fixed external field through a sequence of spin-by-spin flip 
transitions until a new stationary state to be achieved. The main topic of this paper is to 
represent the alternative analytic method that gives a possibility to establish the direct 
connection between the infinite number of metastable states in RFIM and magnetization 
hysteresis loop. The method is based on the calculations of the total number of energy 
minima for the random field Ising Hamiltonian corresponding to different fixed values 
of the external field and specific magnetization per site. The hysteresis looop appears 
as a marginal curve separating m — h plane into the metastability region inside of the 
loop and the region of completely unstable states outside of it. 

2. Local Energy Minima in RFIM 

2.1. Stationary States of Zero Temperature RFIM Hamiltonian 

Consider the Ising spin system (jj = ±1 defined on a regular A^-site lattice and interacting 
with the nearest neighbors, external uniform fild h and irregular Tj (random) fields (here, 
i G 1,...,A^ is a lattice site number). The energy of such a system is known to be 
represented as follows: 

H{s} = -]-J VijCTiCrj- X! + ^» (1) 
ij i 

where, the first term represents the nearest neighbors interaction energy {vij = 1 7^ 
only for nearest neighbors), J is a spin-spin interaction constant and the second one 
reflects the influence of both the external and the random fields; {s} = {ai, o"2, ajq} 
denotes iV-spin configuration vector. The random field values Xj are assumed to be non- 
correlated at the different lattice sites and characterized with the probability distribution 
function P (r) that is the same for all the lattice sites. One can also define the possible 
local minima of RFIM Hamiltonian as a set of Ising spin configurations s satisfying the 
following conditions: 

H{s}<H{s^^} \/iel,...,N (2) 

where, each local minimum s is a spin configuration having lower energy in respect to 
all other configurations s*^*-* that differ from the initial one s by a single spin at i-th 
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site flipped in opposite direction: cTj ^ — aj. In other words, the system of inequalities: 



u,{s} = -[H{s^'^-H{s}] = a,\^T, + h + JJ2vij^jj >0; V2el,...,iV (3) 

define all the possible stationary spin configurations for zero temperature RFIM. It 
should be noted that these inequalities can be satisfied in the unique case when both 

factors in the product (Ji + h + J J2 Vijaj^wiW have the same signs. Therefore, the 

solution of inequalities is completely equivalent to the following equations: 



ai = sign + h + ^ijaj j (4) 

One can interpret these equations as local equilibrium conditions for each site spin in 
the local self-consistent field 

hi = Ti + h + Vijffj (5) 
i 

Accordingly, each site spin is oriented such a way to have the same direction as the local 
field applied. This equation also represents a zero temperature limit of the mean field 
equation that is often used in different studies of the equilibrium properties of RFIM. 
In particular case J = the Eq. (|^) immediately gives the only possible solution: 

(Tj = sign {n + h) (6) 

that corresponds to a single minimum of RFIM in this case. The total magnetization 
per site in this case can be expressed as 

'mo{h) = ^ J2 = J dTP{T)sign{h + r) (7) 

where, denoting the Dirac's function as 6{t) one can define the probability distribution 
function of the random field values: 

m = ^ E ^(^ - (8) 

i 

mo(/i) will therefore represent a single magnetization curve in case of non- interacting 
Ising spin system in a random field. 



2.2. Numerical Simulation of Magnetization Curves 

The most important effect resulting from presence of the infinite number of local 
energy minima in RFIM can be expected for the magnetization properties of Ising 
spin system in zero temperature case when thermal fluctuations do not play any role. 
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In absence of thermal fluctuations each stationary state will still remain the same as 
long as the first external field value is achieved where this state becomes unstable and 
pinning equations Eq.(Q) or Eq.(^) can not be more satisfied. Then, the unstable spin 
configuration is changed at fixed external field through a sequence of single spin flip 
transitions until a new stationary state to be achieved. Therefore, due to great number 
of possible stationary microscopic configuration consistent with the pinning conditions 
the magnetization process can be considered as a bifurcation sequence of transitions 
between the different stationary states as the external field is changed. The necessary 
condition for a single spin transition can be represented as follows: 

(Tj =^ — CTj if only Ui{s} < & Ui{s} = min (9) 

It is denotes that new stable spin configuration at each fixed external field is 
achieved through a sequence spin by spin flip processes running along a single possible 
configuration path providing the maximal energy gain at each elementary stage of 
transition. The configuration path satisfying conditions Eq.(|^) provides the maximal 
rate of the free energy decrease in agreement with usual requirements of non-linear 
thermodynamics and give necessary algorithm for the computer simulation methods 
and modeling of magnetization processes in framework of RFIM. Fig.l shows the typical 
magnetization loop resulting from such a type simulation procedure performed for the 
one-dimensional RFIM. Numerically simulated gaussian random field values on the one- 
dimensional lattice with 1000 sites were used for calculations. 



3. Analytic Accounting of the number of local minima 

3.1. General expressions 

By using the definition from Eq.(^one can propose the following general way for 
evaluation of the total number of local energy minima in RFIM. For this aim, denote 
initially the half of energy difference as: And then, one can easily construct the 
characteristic functional containing information on the all possible local minima in RFIM 
that can be written as follows: 

L{s}=nr]Ms}) (10) 

i 

where, rj (t) is the Heaviside function (?] (t) = 1 , if t > and rj (t) = 0, if t < 0). 
One can easily check that the above introduced characteristic functional takes the unit 
value only if any spin configuration is a local minimum of RFIM Hamiltonian and zero 
in other case. Therefore, taking the sum over all possible 2^ Ising spin configurations 
one can obtain the following general representation for the total number of local energy 
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minima in RFIM: 

z{h)=j:L{s} (11) 

{s} 

Respectively, one can also find the averaged over the random field ensemble number of 
local minima as follows: 

{s} 

Due to statistical independence of the random field values for different lattice sites 
the averaging procedure in Eq.(0) is transformed into the product of one-site averages: 



1 1 



Here, we used the following evident property of Heaviside function: t] (t) + r] (— t) = 1 
and 7] (t) — T] (— t) = sign (t) and also rj (ait) = | (1 + aisign (t)) where, 

mo (t) = J sign {t + r) P (r) dr. (14) 

It will be shown below that mo (t) represents a magnetization curve of RFIM in case 
of non-interacting spins when J = 0. Therefore, the problem of the total local minima 
number estimation for the Random field Ising model can be always reduced to a standard 
statistical-mechanical procedure on calculation of the partition function for a generalized 
Ising model: 



Zih)=Y^exp{V{s}) (15) 

{s} 

with some effective Hamiltonian: 

V^{s} = -N In 2+ ^ In 1^1 + aiUiQ (^h + J J2 ^u^i j j (16) 

It is important that the effective Hamiltonian can be always expressed as a finite order 
polynomial of Ising spin variables: 

V{s} = V'+ Y: V^a,+ Y: V^^a,a,+ ^ Vij,a,a,a, + ... (17) 

i ij ijk 

Polynomial coefficients are translation invariant functions of their arguments. These 
take zero values for coinciding pairs of their arguments and contain only the short range 
nearest neighbor interactions. The maximal order of this polynomial is k + 1 where, 
A; is a number of nearest interacting neighbors of the initial RFIM Hamiltonian. As 
a result, using the analogy with the partition function of equilibrium thermodynamics 
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one can assume the following universal asymptotic behavior for the average number of 
local minima when the lattice sites number becomes infinitely large: 



Z{h) exp(^(/i)iV) (18) 

where /i {h) is the independent on function of the external field. Because the total 
number of local energy minima 1 < Z (h) < 2^, then < n{h) < ln2 and one can 
make a general conclusion that Z {h) must exponentially grow as cxd. Therefore, 

RFIM can have a single minimum if only /i {h) = and infinite number of local minima 
in other case. Such a case is realized for the non-interacting spin system at J = 0. One 
can easily obtain this result by using Eq.(|T3|) at J = 0: 

N 




Z{h)j=o= o 22 a + ^^o{h))\ =1 (19) 



Naturally, in general case both Z (h) and /i (h) are the functions dependent on the 
parameters of Random Field Ising Hamiltonian such as the external field h, spin-spin 
interaction constant J and also on the distribution function of random field P{t). Along 
with the total number of local energy minima Z {h) one can also consider another 
important quantity: Z (M, h) that gives the number of stationary states corresponding 
to a fixed value of total magnetization M =Y^ ai. Similar to Eq.(|T5D this one can be 

i 

found as follows: 



Z{M,h) = Y. exp(V{s}) (20) 

{s;M} 

where, the summation is performed over configurations with a fixed total magnetization 
value. One can also expect the following universal asymptotics in the thermodynamical 
limit N oo: 



Z (M, h) — > exp (/i (m, h) N) (21) 

where, m = M/N is a specific magnetization per site. But unlike to always positive 
function /i(/i) defined above the new function fi (m, h) introduced here and dependent on 
two variables m and h can take the negative values too in the corresponding area of their 



change. As a result of Eq.(21), the number of possible stationary states corresponding to 
different fixed values of specific magnetization m becomes infinitely small approaching 
zero value as A^ — > oo in that region oi m — h plane where /z (m, h) < 0. Therefore, one 
can conclude that all possible infinite set of stationary solutions is always located only 
inside the region /i (m, h) > 0. The marginal curve /x (m, /z) = separating these two 
regions can be therefore associated with a global magnetization loop. Such interpretation 
is completely consistent with the results of dynamical approach based on the computer 
simulation of partial magnetization sub-loops that are always locked inside of main loop 
confirming its marginal character as shown in Fig.l. 
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Due to above mentioned analogy with the equihbrium statistical mechanics a great 
number of approximation methods can be applied to calculation of Z (K) and Z (M, K). 
As an example, the simple perturbation method will be used below to estimate the 
number of local minima in case of the small interaction constant limit . The main reason 
of that approach is to show that the infinite number of local minima effect immediately 
appears as only J 7^ 0, starting from an arbitrary weak spin-spin interaction. 



3.2. Perturbation Theory Expansion 

As known the perturbation theory methods in statistical mechanics providing correct 
behavior in the thermodynamic limit when N ^ 00 are realized as a free energy 
expansion that is proportional to the log value of a partition functional. Similarly, 
the corresponding expansion in respect to In (/i)) = Ni^{h) and ln(z(M,/i)) = 
Nfi (M, h) should be performed in case J 0. Accordingly, one can easily obtain the 
Taylor's series expansion for the effective Hamiltonian defined in Eq.(|l6D as V^{s} = 
Vo{s} + Vi{s} + V2{s} up to the second order on J. Then, using Ising spin algebra rule 
((T^ = cr), a zero order term can be expressed as follows: 

Others first and second order perturbation contributions are: 



^ ^ 2^(1 + a.mo {h)f I Y / 2 ^ 1 + ^«^o {h) \ ^ 



2 {l-ml{h)Yt^ 



1 Tfi { 

^ i - mo (^ftj •^•fc 

As a result, the corresponding expansion of both the In [Z{h)) = N^i{h) and 

In [Z (M, h)) = Nfi (m, h) up to the second order on the interaction constant J will 
look like as follows: 

In (Z{h)) = In (Mh)) + (V,), + {V2), + \ {{V^)^ - (yi)J) (25) 
In {zWJi)) = In (^Z,{M,h)) + {V^),, + {V,),, + l {{v,') - {V^)l,){26) 
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Here, the averaging procedures are defined as: 

E {...)expVo{s} E (...)expK){s} 

/ \ = i!l ■ ( ) = ^ (27) 

{s} {s,M} 

where, {s,M} is a set of spin configurations with a fixed magnetization value M. Due 
to statistical independence of spin variables in the different sites and evident result: 
— ^o{h) that can be easily found from the Eqns. (p2| , pTj) both the {Vi)^ and 



second term from the Eq.(24) contributing into (V2)q take zero values. Taking also into 
account that acorrding to Eq. ([T9|) ZQ{h) = 1 and ln{Zo{h)) = one can obtain the 
following results for the total number of stationary states Z (h): 

1 2 



In [Z{h)) = - J^JT TTT^ Y (^^j) ~ ^^^^ (^))o " ^j^o W) 

^ ^ 2 (1-mUh)) u 



In [Z{h)) = --J^—^ r-7-2 ^ij^ik - mo (h)) (r^ak) + 

2 (l-m^(/i)) ° 

1 2 ('^o (^)) 

7^ ^TTT^ 51 ((^i - "^0 (/i)) ((Tz - mo (/i)) (Tj(Tfe)o (28) 

2 (l-m^(/i)) ij^fc 

The first term in square braces is completely reduced with a part / = i of the second 
sum. One can also check that only the k = ihl = j terms from the remaining part of 
the second sum will give nontrivial contributions. Therefore, 

, 2 

, 2- E 

ml (h)) ij 

\j' H (^))' E KO' = \j" H {h)fkN (29) 

where k is the number of nearest neighbors. Therefore, the total number of energy 
minima becomes infinitely large even in case of an arbitrary weak Ising spin interaction 
and can be expressed in the second order perturbation theory approach as follows: 

Z{h) = exp (/i (h) N) where, {h) = ^kJ^ (m'o (h))^ (30) 
3.3. Multiple minima and magnetization loop 

In this section the relationship between the great number of possible stationary state of 
RFIM and the irreversible magnetization loop typical for ferromagnetic materials will 
be a subject for study and discussion. It will be shown that very important information 
on the magnetization loop can be obtained not only from the direct dynamic simulations 
of the magnetization process but also by accounting the number of possible stationary 
solutions and their distribution on the different values of macroscopic magnetization. 
More exactly, the method is based on the calculation of Z (M, h) and solution of 



characteristic equation for the hysteresis loop: In {^Z{M, h)j = /i (m, h) N = 0. In the 
framework of perturbation theory /i (m, h) can be found directly from the Eq.(^). But 
there is a simplest way to do that by using the evident sum rule: Z (h) =J2 Z (M, h) 

M 

following from definitions of both these quantities. Using this rule and exponential 
representations (see Eqns.(0) and (|21|)) one can obtain asymptotically (as ^ oo ) 
the following results: 



exp (/i (h) N) exp (/i (m, h) N) — exp max (/i 



m,h))N\ (31) 



Or, equivalent relationship: 

^ d _ 

IJ,(h) = n (m (h) , h) where, — — /i (m, h) = at m = m (h) (32) 

dm 

As a result fi {m, h) can be expanded near the maximum m = m (h) as follows: 

1 _ 2 

li{m,h) = —-(3 {h) {m — m{h)) +fi{h); where, P (h) > (33) 

In the case of weak interaction J this representation can be used to find the 
magnetization curve as a solution of the equation fi{m,h) = 0. According to Eq. (pO|) 
/i (/i) = at J = and there is a single degenerated solution m = mo (/i)for a zero order 
equation: 

/io(m,/i) = -^/5o(/i)(m-mo(/i))' = (34) 

So, there is no hysteresis at J = 0. But, as follows from the Eq . (|33|) , it is immediately 
appears as only J ^ and represented by two branches m = m± (h) of hysteresis loop 
close to m = rh {h) and each other, where: 

m± {h) =m{h)± y^2/i {h) /(3{h); where, P (h) > (35) 



Because y /i (h) is the first order on J then only up to the first order calculations for 
m (h) = tuq (h) +mi (h) and zero order for /3 (h) = Pq (h) should be performed. In 
particular, one can easily show that m {h) is exactly equal to a single-site average of 
Ising spin value (ct;) and expanded into the perturbation series as: 

m{h) = {(Ji) = {cri)^ + {{(7i - {<7i)q)Vi (s))^ = mo {h)+kJm'Q {h)mo {h){36) 

Po {h) can be found from the direct calculation of zero order term /zo {rn, h) = 
In (^Zo(M, /i)^ according to Eq.(^) at J = and its expansion near the maximum 
at m = mo (/i). Therefore, Po{h) = (— <9^/c^"^^) /^o (^5 ^)|m=mo(/i)) • Results of such 
calculations are represented below: 



Zo{M,h)=Y^ expiVois}) (37) 

{s;M} 
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/io (m, h) = j^ln (Zo(M, h) 



-[(1 + m) In 



1 + mo (/i)' 
1 + m 



+ (1 — m) In 



'1 - mo (/;.)" 



m 



(38) 



Po {h) = (-8^1 dm^) fio (m, h) = ^ _ (39) 

The final equation representing the irreversible magnetization loop for a zero 
temperature RFIM in the weak interaction limit will look like as follows: 



m± (h) = (l + Urn'o (/i)) mo (h) ± Jm^ (h) ^k{l-ml{h)) (40) 

A set of hysteresis loops calculated according to Eq.(|y) for different values of the 
interaction parameter J is indicated in Fig.2. 

Therefore, the present approach makes it possible to perform detailed anlytic 
studies concerning physical nature of hysteresis by using the Random Field Ising based 
microscopic consideration. 
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Figure cuptions 

Fig.l Magnetization loop and partial sub-loops numerically simulated for the one- 
dimentional RFIM. 

Fig. 2 Marginal hysteresis loops separating metastable band unstable states in 
RFIM calculated for the different values of interaction constant in the framework of 
perturbation method. 
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